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A COMPLETE SURVEY OF CASE A BINARY EVOLUTION WITH COMPARISON TO OBSERVED 

ALGOL-TYPE SYSTEMS 
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ABSTRACT 

We undertake a comparison of observed Algol-type binaries with a library of computed Case A binary evolution 
tracks. The library consists of 5500 binary tracks with various values of initial primary mass Mio, mass ratio qa, 
and period Pq, designed to sample the phase-space of Case A binaries in the range -0.10 < logMio < 1 .7. Each 
binary is evolved using a standard code with the assumption that both total mass and orbital angular momentum 
are conserved. This code follows the evolution of both stars until the point where contact or reverse mass transfer 
occurs. The resulting binary tracks show a rich variety of behavior which we sort into several subclasses of Case 
A and Case B. We present the results of this classification, the final mass ratio and the fraction of time spent in 
Roche Lobe overflow for each binary system. The conservative assumption under which we created this library 
is expected to hold for a broad range of binaries, where both components have spectra in the range GO to Bl and 
luminosity class III - V. We gather a list of relatively well-determined observed hot Algol-type binaries meeting 
this criterion, as well as a list of cooler Algol-type binaries where we expect significant dynamo-driven mass 
loss and angular momentum loss. We fit each observed binary to our library of tracks using a x^-minimizing 
procedure. We find that the hot Algols display overall acceptable x^, confirming the conservative assumption, 
while the cool Algols show much less acceptable x^ suggesting the need for more free parameters, such as mass 
and angular momentum loss. 

Subject headings: stellar evolution 



1. INTRODUCTION 

Many binary stars are observed to be undergoing Roche-lobe 
overflow (RLOF), which is recognised as being a natural re- 
sponse to the fact that, for a binary of given separation, there 
is a critical maximum radius, the Roche-lobe radius, that a star 
cannot exceed without losing mass to its companion. There 
are many sub-types of stars undergoing RLOF, but we concen- 
trate here on those which, like the prototype Algol, consist of 
(i) a lobe-filling, mass-losing star that is substantially above the 
main sequence, and (ii) a component which underfills its Roche 
lobe, and is usually nearer to, though still larger than, the main 
sequence. We concentrate on those (Case A) with short initial 
periods, the lower and upper period depending on the primary 
mass. 

It is not difficult to evolve theoretically pairs of stars with 
a given initial primary mass Mio, initial mass ratio ^o and ini- 
tial orbital period Pq, and follow them into, and beyond, the 
stage of RLOF. However, such evolution is certainly affected 
by assumptions regarding both mass loss and angular momen- 
tum loss from the system as a whole. As a zero-order model 
it is commonly supposed that both total mass and orbital angu- 
lar momentum are conserved, and we have computed conserva- 
tive evolution for a large number of binary initial parameters: 
37 X 10 X 15 models with various Mio, ^o and Pq. Most of the 
periods considered are appropriate to Case A, but some corre- 
spond to Case B. 

There is plenty of evidence, both direct and indirect, that 
mass loss and/or angular momentum loss takes place in at least 
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some systems. If mass escapes from the system as stellar wind, 
then it will also carry angular momentum away. Mass loss is 
observed fairly directly both in cool stars, where it appears to 
be driven by dynamo activity in their convective envelopes, and 
in hot stars, where radiation pressure in spectral lines may be 
the main driving force. Mass loss is also clearly evident in many 
stars of supergiant luminosity, across the whole range of spec- 
tral types; but we do not consider supergiants here. However 
there is a broad range of spectra, from about GO to perhaps Bl 
and luminosity class III - V, where there is rather little evidence 
of significant mass loss, and where the conservative assumption 
may therefore be reasonable. We test this by comparing a selec- 
tion of observed 'hot Algols' (having both spectra in this range) 
with theoretical conservative models, using a x^ test. We find 
a reasonable agreement, especially if we exclude one system 
which is near the extreme of this temperature range. Comparing 
the same conservative models against some observed 'cool Al- 
gols' we find, as we expect, that the agreement is much poorer. 
We have used a massively-parallel array, the Compaq Tera- 
cluster 2000 at LLNL, to evolve our data cube of models. This 
data cube covers the following ranges of initial primary mass 
Mio (in solar units), initial mass ratio, defined by 

Mio , 

and initial period Po: 

logMio = -0.10,-0.05, 

log^o = 0.05,0.10, 

log(Po/i'zAMs) = 0.05,0.1,. 
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,1.7, 


(2a) 


,0.5, 
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0.75. 
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Here Pzamsi ^ function of Mio, is the period at which the ini- 
tially more massive component would just fill its Roche lobe on 
the zero-age main sequence. We used the approximation 



fZAMS ■ 



0.19Mio + 0.47Mfo* 
1 + 1.18M2„ 



2.33 



(3) 



These initial periods cover Case A and a small part of Case B. 
We constructed such a 'data cube' with each of six metallicities 
(Z = 0.03 , 0.02, 0.01 , 0.004, 0.001 , .0003), and also, for Z = 0.02 
only, with three different assumptions about mass loss/angular 
momentum loss (in addition to the conservative assumption). 
We present here only the conservative, Z = 0.02 data cube. 

In §2, we discuss the numerical modelling and the physical 
assumptions that go into our data cube, and in §3 we discuss the 
results. We attempt to classify the results into a small number of 
sub-categories of Case A (and some analogues in Case B), de- 
pending for instance on whether the two components come into 
contact rapidly, slowly, or not at all after the start of RLOF, and 
(in the last case) on whether or not primary reaches a supernova 
before the secondary swells up enough to reach reverse RLOF. 
In §4 we discuss our attempts to fit several observed semide- 
tached systems (Algols) with the theoretical models. We give 
our conclusions in §5. 

We emphasise here that even if a particular Algol can be 
reasonably fitted by a conservative model, this does not prove 
that the evolution was conservative. Some models of non- 
conservation might lead to the same current parameters, start- 
ing from different initial conditions. Even if we had a mass 
loss/angular momentum loss model with no free parameters 
in it, we might still have ambiguity, partly because there 
are only six independent observational parameters (current 
P, M\ , M2, R2, T\, T2) to be fitted by four theoretical parame- 
ters (age, and Pq, Mio, Io), and partly because our data cube is 
still quite coarse even with 5550 models in it. 

We also emphasise that throughout this paper we use suffixes 
1 and 2 consistently to refer to the components with the greater 
and smaller initial mass respectively. This may seem unfor- 
tunate since observers normally call the currently hotter (and 
normally more massive) component the 'primary', at least in 
Algol systems. This component is the descendant of the origi- 
nally less massive star We do not think it would be helpful to 
interchange the suffices at the points in evolution where the or- 
dering of the temperatures changes. However to avoid the most 
obvious possibility of confusion we do not use the terms 'pri- 
mary' and 'secondary': instead we refer to the components as 
*1 (pronounced 'star one') and *2, and keep these designations 
throughout their entire evolution. The mass ratio q as defined 
in Equation 1 always starts off with qQ> 1. After some RLOF 
it is commonly < 1 . 

2. THE THEORETICAL DATA CUBE 

We used the stellar evolution code most recently described by 
Pols et al. (1995), based on the code of Eggleton (1971, 1972); 
Eggleton, Faulkner & Flannery (1973). This code is fully im- 
plicit in the composition equations as well as in the structure 
and the mesh-spacing equations. The implicit adaptive mesh is 
particularly useful for mass-transfer situations. In fact, it means 
that in a first approximation we do not have to do anything to 
the code to account for mass transfer, except replace a bound- 
ary condition M{t) = constant by a condition which gives the 
mass-loss rate M as a function of stellar radius R and Roche- 
lobe radius i?L (Tout & Eggleton 1988). 



We will not repeat here a description of the physical input 
(Pols et al. 1995). We have however included a simplistic 
model of convective overshooting (Schroder, Pols & Eggleton 
1997; Pols et al. 1997), based on a comparison of theoretical 
and observed non-interacting binaries. Other assumptions in 
the code are standard, and include the following: (a) the con- 
vective mixing of composition is treated as a diffusion equa- 
tion, with diffusion coefficient a function of Vr - Va (Eggleton 
1972), and (b) because the mesh is fully adaptive, i.e. non- 
Lagrangian, an upstream advection term is needed in all time 
derivatives (Eggleton 1971). The former ensures that any con- 
vection zones satisfy the K. Schwarzschild convection crite- 
rion (Vr ~ Va) and simultaneously that any semiconvection 
zones that may arise satisfy the M. Schwarzschild condition 
(Schwarzschild & Harm 1958) and are dealt with automati- 
cally, without extra code; the latter ensures that any evolution- 
ary stage involving thin burning shells is computed very effi- 
ciently. 

Regarding situations specific to binaries, we make the fol- 
lowing assumptions: 

(i) The star is still treated as spherically symmetric, the 
radial coordinate r being the volume-radius of an 
equipotential surface. The gravity at effective radius r 
is reduced by a factor dependent on angular velocity: 



Gm 



3Gm 



(4) 



(5) 



where m is the mass within an equipotential of 
volume-radius r and il is the angular velocity of the 
star, assumed corotating with the binary. 

(ii) Mass transfer from a star that overfills its Roche lobe is 
treated as spherically symmetric, and governed by the 
boundary condition 

M =-C{log(R/RL?}, R>Rl 
= 0, R<Rl 

where C = 5OOM0/yr Thus a transfer rate of 
5 X lO~^M0/yr corresponds to an overfill of 0.1%. We 
only do this for *1. There inevitably comes a point 
in evolution when *2 fills its own Roche lobe, but 
this is usually either (a) while * 1 already fills its own 
lobe, so that the binary comes into contact - for the 
present, we stop evolution at this point; or (b) when * 1 
has evolved to a late and relatively compact state of 
low mass, and *2 has grown to a very large radius. In 
the latter case the mass ratio is very small (q < 0.2). 
Therefore the mass transfer can be expected to be rapid, 
and unstable on a short (hydrodynamical) timescale. 
We expect common-envelope evolution beyond this 
point (Paczyriski 1976), and so we stop evolution at this 
point, also. 

(iii) It is assumed that the matter which leaves * 1 is accreted 
in a spherically symmetric manner at the surface of 
*2, with entropy and temperature equal to the surface 
values of *2. Thus no model is incorporated for the 
temperature/entropy budget of the material during 
transfer. This may seem potentially serious, but when 
most of the mass is transferred on a nuclear timescale it 
should not be important. 
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(iv) The composition of accreting material on *2 is assumed 
to be the same as that of material already just below 
the surface of *2, rather than (as it should be) of the 
material leaving *l. This is just done for convenience, 
and is only significantly in error at a fairly late stage 
in mass transfer The observed Algols that we make 
comparison with are probably not at such late stages. 

(v) On a somewhat technical level, the implementation of 
Equations 4 and 5 numerically within the framework 
of a fully implicit and adaptive code means that it is 
desirable, though possibly not essential, to introduce an 
extra equation into the usual set of difference equations 
- for the structure, composition and mesh-distribution 
variables - that are solved for by Newton-Raphson 
iteration. This is because Equation (4), while depending 
primarily on the local variables r and m(r), also 
depends on the surface mass M(t), through il(0- ^ 
depends not only on the orbital angular momentum 
(which, being assumed constant in a conservative 
model, is no problem) but also on the the masses of the 
two components via Newtonian gravitation. Because 
of Equation (5), M{t) is not known a priori, but only 
after the iteration is finished. We found it convenient to 
add M'(t,r) as a new but somewhat artificial variable 
satisfying the trivial equation 

f^=0, (6) 

or 

with the equally trivial boundary condition that at the 
surface 

M'(t,r) = M(t). (7) 

Although this modification is barely necessary for the 
conservative models, it is rather more important for 
non-conservative models, where Equation 5 may have 
an extra term, attributable to stellar wind, and where 
the angular momentum is no longer constant. 

With the above assumptions and modifications the code 
works reasonably satisfactorily in an automatic way. We set up 
a grid of starting models with Mio, log qo and Pq given by Equa- 
tion 2. For most masses in the range I.SMq < M < 16Mq we 
found that we obtained Case A evolution with 1 < Pq/^zams ^ 
4, and Case B for > 4, but the critical value for Case B de- 
creased rapidly below ^ 1.5Mq, and increased slowly above 
- 16Mq. 

Given Mio, qo and Po, we started by evolving *l until one or 
other of the following conditions occurred: 

(a) 2000 timesteps were taken 

(b) carbon-burning luminosity exceeded ILq, indicating 
that a supernova explosion was imminent 

(c) the age exceeded 20 Gyr 

(d) the code failed to converge, or 

(e) the stellar radius exceeded the Roche-lobe radius by 
more than 10%. 

For * 1 , (e) indicated that hydrodynamically unstable RLOF 
was taking place, usually due to a large initial mass ratio 
(^0 ^ 3), or to a deep convective envelope on the loser. 



We then ran *2, giving it a rate of mass gain that was the 
negative of the stored mass loss rate of * 1 . This run was also 
terminated at the first point when one of the first four conditions 
above occurred, but it could also terminate itself if 

(f) the age of *2 went beyond the age at which * 1 
terminated, or 

(g) the radius of *2 reached its Roche-lobe radius. The 
latter normally meant either that the system had evolved 
into contact, *2 filling its lobe while * 1 still was, or else 
that it had evolved into a reverse RLOF situation, with 
q < 0.2, that would presumably lead to mass transfer 
on a hydrodynamical timescale, probably implying 
common-envelope evolution. In either case, the 
implicit assumption that * 1 's evolution is independent 
of whatever happens to *2 breaks down, and so we 
consider here only the evolution that takes place prior 
to the point where *2 filled its Roche lobe. 

Convergence failure - (d) above - was not very common, 
though more common than we would have wished. For *l 
it was usually because either (i) Equation 3 apparently gives 
slightly too small a value, for some ranges of Mio, so that at the 
lowest value of Pq for those masses * 1 already filled its Roche 
lobe while still making a rapid adjustment from the approxi- 
mate ZAMS from which it started; or (ii) for the most mas- 
sive stars, > 25Mq, a breakdown often occurred when *1 ap- 
proached a sloping line across the HRD, starting just before the 
terminal main sequence at ^ 50Mq (our highest initial value) 
and reaching to the red supergiant region at ^ 25Mq. It may 
not be just coincidence that this is also approximately the ob- 
servational 'Humphreys-Davidson Limit', which appears to be 
an upper limit for stars in the HRD. Stars close to this limit 
are typically P Cyg stars, Hubble-Sandage variables, or Lumi- 
nous Blue Variables (LBVs). Such stars have internal lumi- 
nosities that are close to or even above the Eddington limit in 
zones where the opacity has a local maximum. Thus it may be 
that the numerical convergence difficulties have their origin in 
the physical difficulty of maintaining hydrostatic equilibrium in 
such stars. 

For *2 convergence failure occurs because: 

(i) In binaries with extreme initial mass ratios q^ > 4, 
*2 often fails to converge while gaining mass at the 
thermal timescale of *L This may occur because the 
large mass ratio means that the thermal timescale of *1 
is closer to the dynamic timescale of *2 

(ii) In our lowest mass binaries, M20 ;$ O.8OM0, the mass 
gaining star (*2) 'ignores' the fact that *1 is losing mass 
for a handful of timesteps and maintains a constant 
mass. Then *2 attempts to gain all the mass that *1 has 
lost in of order ten timesteps in a single timestep and 
fails to converge 

(iii) When the mass gaining star has a mass ^ 1.5Mq. 
These stars are in the transition region between lower 
MS stars with convective envelopes and upper MS stars 
with convective cores, and they possess very shallow 
surface convection zones which may only be a few 
mesh points wide. We suspect that this barely resolved 
surface convection zone contributes to their numerical 
instability; or, finally. 



(iv) we also see a theoretical 'Humphreys-Davidson Limit' 
in *2 at high mass. 

3. CLASSIFICATION OF TYPES OF EVOLUTION 

We define here the six major subtypes of Case A evolution 
identified by Eggleton (2000), cases AD, AR, AS, AE, AL, AN. 
In addition, we define two rather more rare cases, AG and AB. 
Three of these subtypes (AD, AR, AS) lead to contact while 
both components are on the main sequence (MS). Two cases 
(AE, AG) reach contact with one or both components evolved 
past the terminal MS. After the initial episode of mass transfer 
from *1 to *2, the remaining three cases experience a period 
of separation followed either by reverse mass transfer at very 
small q (AB, AL) or the supernova of * 1 (AN). Specifically, the 
six cases, are: 

• AD - dynamic RLOF: this occurs in binaries with 
large ^o and in binaries where the mass losing star (*1) 
has a deep convective envelope. Once RLOF begins, 
mass transfer quickly accelerates to the dynamic 
timescale of *1, fdyn, which we assume to be less than 
a tenth of the thermal timescale, Ikh- The thermal (or 
Kelvin-Helmholtz) timescale is determined in the code 
as the integrated total energy, thermal plus gravitational, 
divided by the total luminosity at the surface. Thus the 
mass transfer is determined to be dynamic when 



M> 10.0 X 



M 

fKH 



(8) 



The calculation is terminated by (e) above, but seems 
likely to lead either to contact or to a common-envelope 
situation, and probably then to a complete merger of 
the two components. We illustrate the behaviour in the 
HRD and the mass transfer rate of case AD in Figure 1 . 

• AR - rapid evolution to contact: this occurs in 
binaries with moderate to large ^o- In these cases *2 
expands so rapidly in response to the onset of ^I's 
thermal-timescale RLOF that it fills its own Roche 
lobe before much mass is transferred. We define the 
mass transfer rate to be thermal when the magnitude 
of the thermal luminosity, | Ltherm |, reaches 2% of the 
nuclear burning luminosity, Lnuc- This probably leads 
to a contact binary of the W UMa type, although it 
can happen as easily for massive stars (provided qo 
is suitably large) as for the lower masses of typical 
W UMa systems. Case AR behavior is illustrated in 
Figure 2. 

In some binary runs, these two cases are difficult to dis- 
tinguish. While evolution of *1 will proceed through several 
timesteps of dynamic timescale mass transfer before being ter- 
minated by (e), the calculation of *2 is often unable to converge 
while gaining mass at this rate. The calculations of case AR and 
AD binaries at very large qo, therefore, often terminate before 
contact is reached and we must guess the rate mass transfer 
achieves before contact occurs. 

To do so we extrapolate the function log^(f) to the time 
^contact at which the radius of *2 has expanded to fill its RL and 
log^(f) = 0. We then examine the mass loss history of *1 
(whose calculation has proceeded further in time than that of 



*2) and determine whether the mass transfer rate reaches the 
thermal or dynamic timescale at f < fcontact- Unfortunately, the 
function log^(f) can be both non-linear and slightly noisy, 
and so fcontact, along with the maximum mass transfer rate, can 
depend rather sensitively on the exact point in time at which *2 
fails. 

• AS - slow evolution to contact: this occurs in binaries 
with small qo and small Pq. These binaries experience a 
short burst of thermal timescale mass transfer, followed 
by a long phase of nuclear timescale mass transfer, 
during which much mass is exchanged. The two stars 
come into contact slowly, but reach contact before 
either star has left the MS. The large amount of mass 
transfer leads to a final mass ratio substantially below 
unity (typically q ^ 0.4-0.6), and with both stars 
substantially larger than their ZAMS radii. Case AS 
behavior is illustrated in Figure 3. We note that while 
*2 always remains near the main-sequence band, * 1 
evolves to substantially cooler temperatures. This is a 
common configuration in observed Algol systems. 

• AE - early overtaking: this occurs in binaries with 
small ^0 and moderate Pq- It occurs only in binaries 
with initial masses 2Mq <Mi < IOM0. The mass 
transfer in this case is very similar to case AS. In 
case AE, however, *2 gains so much mass that its 
evolution is accelerated to the extent that *2 reaches 
the Hertzsprung Gap, HG, while *1 is still on the MS; 
the evolution of the initially less massive star, *2, has 
overtaken that of * L We define the overtaking as early 
because it occurs with *1 still on the MS. Most case AE 
binaries reach contact shortly thereafter. However, in a 
few cases * 1 shrinks very slightly inside its RL at the 
end of the calculation and the run ends with the RLOF 
of *2. In these cases, * 1 has very nearly exhausted 
hydrogen and it is likely that it will soon swell once 
again to fill its RL and contact will again occur. Case 
AE behavior is illustrated in Figure 4. 

In most cases where contact is avoided while *1 is on the 
MS, *1 loses so much mass that it eventually shrinks inside its 
RL leaving only a compact core. A period of separation en- 
sues which may then be followed by further RLOF of * 1 or *2. 
These are the cases AL, AB and AN, described in more detail 
below. However, in our lower mass binaries {M\o < 1 .6Mq) we 
see a few cases where contact is avoided while * 1 is on the MS, 
but reached later on. 

• AG - contact on giant branch: this occurs for 

Mio < 1 .6Mq, and Po larger then those of AS/AE, but 
smaller then AL/AN. Contact is avoided while *1 is on 
the MS, but occurs when *1 reaches the giant branch, 
GB. At time of contact *2 is in the HG or on the GB 
as well. A typical example of case AG is shown in 
Figures. 

Cases AL, AN are distinguished by whether or not *1 super- 
novas before *2 reaches RLOF. In practice, we assume a super- 
novae explosion to be iminent when *l begins burning carbon. 

• AL - late overtaking: this occurs in binaries with 
Mio < 13Mq and moderate to large Pq. In these 
binaries, *2 reaches RLOF before * 1 begins burning 
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carbon. In many of the lower-mass AL cases, * 1 has 
become a low mass remnant (WD or NS) which will 
never supernovae unless the (uncomputed) reverse 
mass transfer results in significant mass gain for * 1 . 
The evolution of *2 has overtaken the evolution of *1 
in the sense that the initially more massive star is now 
shrunk inside its RL while the initially less massive star 
is undergoing RLOF. The overtaking is late because 
it occurs with *1 past the MS. Case AL behavior is 
illustrated in Figure 6. 

• AN - no overtaking: this occurs in higher mass binaries 
with moderate to large Po- In these binaries *1 reaches 
carbon burning, indicating an imminent supernova, 
before *2 has reached RLOF. Case AN behavior is 
illustrated in Figure 8. 

As discussed in §2 the evolution of our most massive stars, 
M\o > 25Mq often breaks down. This leads us again to the situ- 
ation where we must make a best guess as to what happens after 
the run stops. To distinguish case AL from AN we must deter- 
mine whether *2 reaches RLOF before * 1 ignites carbon. This 
procedure is somewhat uncertain and leads to the great major- 
ity of our unclassified runs at very high mass. We also suspect 
that several of the highest mass runs, Mio > AOMq, which were 
classified as AL are uncertain, and may more probably be Case 
AN. (See Figure 9.) 

We note that the definitions of case AL and AN given here 
do not correspond exactly to those of Eggleton (2000). In this 
previous work, case AN also included those binaries where * 1 
had become a WD or NS before *2 filled its RL. In this work, 
those binaries are included in case AL. 

Pols (1994) noted that occasionally, in what we call case AL 
here, *2 could get to carbon ignition before *1, and thus be 
the first component to explode as a supernova. Pols (1994) 
modeled, in a simple way, the effect of the presumed common- 
envelope phase in ejecting *2's envelope, and then continued 
the evolution of the core. Although it would always have started 
He burning later than *1, it might be sufficiently more massive 
to overtake and ignite carbon first. However, in our work we did 
not attempt to model the common-envelope phase at all, and so 
we cannot be definitive about this possibility. 

In addition, we include one more class: the classic Case AB. 
In our context this is a subclass of case AL, where *1, after 
becoming a compact helium core with a mass of ^ 1 -2Mq, 
expands again and experiences a further period of RLOF. 

• AB - this occurs in binaries with 6Mq < M\q < 1 IMq, 
at small mass ratios and in a narrow range of periods 
between cases AL and AN. During the second burst 
of mass transfer, *1 ignites helium. It shrinks inside 
its RL for awhile, becoming a compact helium star. 

It then expands again and experiences a third period 
of mass transfer. Although these binaries often fail to 
converge at some point during this third period of mass 
transfer, we suspect that it is followed by a period of 
separation and then reverse mass transfer, making this a 
subclass of AL rather then AN. An example of case AB 
evolution is shown in Figure 7. 

In Table 1 we summarize the seven major subcases (exclud- 
ing case AB), providing the defining equations as well the evo- 
lutionary state and geometrical configuration of the binary com- 
ponents at the end of the calculation. In this table we denote 



the main sequence as M, the Hertzsprung Gap as H, the giant 
branch as G, and low and high mass remnants as R and C, re- 
spectively. In addition, we define the time of first RLOF, fRLOF, 
the approximate MS lifetime of a single star, fMS, the time at 
which the star enters the Hertzprung Gap, tu and the time at 
which carbon is ignited fcbum- We emphasize that we execute 
the classification of each binary in our library automatically, 
and while the various clauses we define work for the great ma- 
jority of systems, we inevitably make a few misclassifications. 

As mentioned above, figures 1-8 illustrate the behaviour in 
the HRD of the subtypes of Case A. We also show the mass 
transfer rate for times when M > 0. Figure 9 shows which 
elements of our data cube reached which outcome. Some of 
the systems of longer Pq are Case B rather than Case A. These 
are usually analogous to either AD, AR, AL or AN. Case BD 
is effectively the classical Late Case B, where *1 reaches the 
giant branch and acquires a deep convective envelope before 
RLOF begins; however it can also be an extreme initial mass 
ratio rather than a convective envelope which triggers dynamic 
mass transfer. Case B systems, or at least those which we have 
computed here, normally have fewer options than Case A be- 
cause it is difficult for *2 to catch up with * 1 when * 1 has al- 
ready reached the terminal main sequence before RLOF. How- 
ever, as emphasised by De Greve & Packet (1990), it is possible 
for early Case B systems to show what we call here Case BL 
for late overtaking, with *2 evolving to fill its own Roche lobe 
while * 1 has shrunk inside its own. This kind of behaviour is 
particularly prevalent in the mass range Miq < 8Mq. 

In Figures 10 and 11 we show, also in colour-coded form, the 
following two properties of systems in our data cube: (i) the fi- 
nal mass ratio of each system and (ii) the fraction of time spent 
as a semidetached system. We define the state of the binary 
to be final when a) contact is reached, b) reverse mass transfer 
begins, c) *1 has ignited carbon, or d) the binary is detached 
and we believe reverse mass transfer to be imminent (ie. the 
function log^(f) will soon reach zero). The time spent as a 
semi-detached system has implications for the frequency of Al- 
gols in the field. We find that for a given primary mass the 
longest-lived Algols originate from systems where mass trans- 
fer begins near the transition to the HG (Late Case A to Early 
Case B or cases AL, AN, BL, BN) and with small to moderate 
initial mass ratios. 

4. COMPARISON WITH OBSERVED SYSTEMS 

Many observed binaries are semidetached (Algols), and one 
might hope that they could be matched by some of the above 
theoretical models during their stage of RLOF. However it has 
been clear for some time (Refsdal, Roth & Weigert 1974) that at 
least some systems (specifically, AS Eri) have such low angular 
momentum that they could hardly have started as detached sys- 
tems of two zero-age main sequence stars of comparable mass. 
Furthermore there are some Algols of such low total mass (e.g. 
R CMa) that they also could hardly have started in such a con- 
figuration. 

In an important paper, Maxted & Hilditch (1996) identified 
9 Algol systems for which they thought the observational data 
was of an unusually high quality. They compared these with 
models computed by De Greve (1993). The comparison was 
not at all satisfactory, the theoretical models having luminosi- 
ties at least 20 times greater than the observed models. They 
also had substantially longer periods. These discrepancies ap- 
pear to be due to the following two features: 



(i) the theoretical models were all Case B 



(ii) they were non-conservative, the assumption being 
made that 50% of the mass lost by * 1 escaped to 
infinity, and 50% was accreted by *2. The escaping 
mass was assumed to remove the same specific angular 
momentum as resided in the orbit of *1. 

We feel that although the kind of non-conservation modeled 
by De Greve (1993) may perhaps be appropriate for massive 
stars (O, and even early B), where radiation pressure may be 
an important agent in mass loss, it is not appropriate for mid- 
main sequence stars where, at least in single stars, very little 
mass loss is normally observed. At the other end of the main 
sequence, stellar winds are rather commonly observed, partic- 
ularly in rapidly rotating G/K/M dwarfs (and even more so in 
giants). These winds probably do not carry off much mass (al- 
though see later), but they may be rich in angular momentum 
because of magnetic linkage to the parent star. We therefore 
think that conservative models may be reasonable for systems 
which are in the middle of the main sequence initially (say Bl 
to GO), and where the loser has not yet evolved to the red-giant 
region at spectra type ~G or later. Following Popper (1980) we 
refer to these systems as 'hot Algols' . Unfortunately rather few 
of the Maxted & Hilditch (1996) selection qualify as hot Algols 
in this sense, although two (U CrB and AF Gem) are on the bor- 
der, with the cooler component having spectral type '--^GO. We 
have therefore included a few more from the literature. Our 
selection of hot Algols is listed in Table 2, with references. 

The observed parameters which we attempt to fit with our 
theoretical models are the six independent quantities logP, 
logMi, logq, log/?2, logTi and logTa. Ri is not independent 
of these, since it is obtained from the assumption that * 1 fills 
its Roche lobe, whose radius is determined by the first 3 pa- 
rameters. Li and L2 are similarly not independent of these 6 
parameters. Our theoretical models have four independent pa- 
rameters, logPo, logMio, log^o and age. 

For each system in Tables 2 and 3 we give three rows. The 
first gives the observational data from the literature, and the 
next the theoretical values from our data cube which minimize 
X^. The second row also includes the best-fit age, in units of 
Myr. The third row gives the zero-age values for the system 
which we infer from our best fit. We use mass-ratio q because 
this is usually obtained more directly from the observational 
data, whether spectroscopic or photometric, than either Mi or 
M2. We list observational errors (when available) in the first 
row for all quantities, but we list total errors (described below) 
in the second row only for those quantities that we actually fit. 

In fitting observed stars to theoretical models, a x^ test seems 
appropriate. However, we have to modify the standard test in 
order to incorporate the fact that our theoretical models have 
an intrinsic 'graininess' because they have not been computed 
for a continuous range of input parameters, but only at the grid- 
points in our data cube. We therefore use a total error, a, which 
is the sum in quadrature of the observational error, cTobs, and a 
'theoretical error', crth, representing the intrinsic graininess. For 
logP, logMi and log ^ we take crth = 0.05, the initial spacing of 
our grid. For log/? and \ogT we take the graininess to be the 
difference in these parameters between adjacent ZAMS models 
from the grid, centered on the mass of the observed binary. For 
example, for an observed star of mass logM = 1 .02 we take the 



theoretical error in the radius to be 

ath,R(logM= 1.02) = /?zAMs(logM= 1.05)-/?zAMs(logM= 1.00) . 

(9) 
We can then look in our data cube for the minimum value of 

(obs-th)2 



X 



=E 



'^L+< 



(10) 



We find that the best fit point picked by minimizing this x^ is 
insensitive to the exact definition of <t^. However, the mag- 
nitude of x^j^ depends directly on a, so we have attempted a 
reasonable definition. In Figure 12 we present the residuals to 
the fit for all Algols from both Tables 2 and 3. 

The hot Algols of Table 2 have a mean x^ of ~ 3. Since there 
are 2 degrees of freedom (6 observed parameters less 4 theoret- 
ical parameters), this value is rather more, but not enormously 
more, than is expected for a normal distribution of errors. The 
number of systems which we use is too small to provide a re- 
ally convincing confirmation or refutation. The worst case, AF 
Gem, is very close to the lower temperature limit, where we 
suppose a priori that conservation might break down. If we 
reject AF Gem, we have a mean x^ of just 2. 

After AF Gem, the next worse cases are DM Per and A Tau. 
Interestingly, both of these systems possess a close third body 
- extraordinarily close in the case of A Tau. The latter sys- 
tem can be seen to be problematic even without a detailed at- 
tempt at fitting. The angular momentum of this system is seen 
to be quite low compared with a system of comparable-mass 
stars at the same total mass, so that something like Case AS 
is to be expected. But Case AS normally evolves into contact 
at a mass ratio which is moderately small, roughly > 0.4 (Fig- 
ure 10), whereas A Tau has quite a small present mass ratio 
of 0.27. This suggests that A Tau has lost some angular mo- 
mentum, necessarily during its slow, nuclear-timescale, RLOF 
rather than the comparable interval of detached evolution be- 
fore RLOF. DM Per's problem is similar, though not so obvious 
without a detailed attempt at fitting. 

But for A Tau and DM Per, unlike most other hot Algols, 
there does in fact exist a mechanism that should do just that. 
The third star in the A Tau system (Fekel & Tomkin 1982) is 
in such a close orbit (33d) that it must influence the orbit of 
the eclipsing pair to a small but significant extent, making its 
eccentricity fluctuate by ^ 0.7% on a timescale of days (Kise- 
leva, Eggleton & Mikkola 1998). Tidal friction will tend to 
oppose this, but can only do so by draining energy and angu- 
lar momentum from the short-period orbit. Conservation laws 
require the angular momentum lost by the inner orbit to go to 
the outer orbit, but the energy loss leads to a net secular evo- 
lution, the inner orbit shrinking while the outer widens. This 
process was probably negligible in the pre-RLOF state because 
the orbit would have been substantially smaller than at present, 
at least if ^0 were not unusually large. But it can now be sig- 
nificant as the stars are larger and the inner orbit wider Tidal 
friction should be capable of setting up a transient equilibrium 
between nuclear evolution, leading to expansion of the inner or- 
bit, and tidal friction, leading to contraction (Kiseleva, Eggleton 
& Mikkola 1998). 

DM Per also has a third body in orbit. The outer orbit is 
longer (^ lOOd) while the inner orbit is shorter, and so the 
process might be thought less likely to be significant. But on 
the other hand the third body is relatively much more massive, 
which may compensate to some extent. 

When we turn to a selection of cooler Algols (Table 3) we 
find significantly larger x^ for many systems. This, we believe. 



Case A Binaries 



is consistent with the view that they are less conservative, cer- 
tainly of angular momentum (which is fairly readily removed 
by magnetic braking on something like a nuclear timescale), 
and perhaps also of mass. We do not normally think of stellar 
winds from cool dwarfs and subgiants as being strong enough 
to remove significant mass, and yet certain active (RS CVn) 
binaries show evidence to the contrary. Both Z Her (Popper 
1988) and to a lesser extent RW UMa (Popper 1980; Scaltriti et 
al. 1993) exhibit the phenomenon that the cooler, presumably 
more evolved, subgiant is the less massive star, despite the fact 
that it does not fill its Roche lobe. This suggests that mass loss 
by wind from the cooler star is already on the nuclear timescale 
of the star 

V1379 Aql (Jeffery & Simon 1997) is an example of a 'post- 
Algol' binary: the low-mass SDB component is presumably the 
remains of * 1 after it has retreated within its Roche lobe, and 
*2 has already evolved to the giant branch. We include it in our 
list of cool Algols as we believe the components to have been 
relatively cool during its Algol phase. The cool ZAMS temper- 
atures which we derive support this assumption. Even without 
detailed fitting, it is clear that the system is problematic. For a 
period as short as 21d the SDB mass is rather large - such a core 
would seem to imply a period of 50 - lOOd. More intriguingly, 
the orbit is very significantly eccentric: e = 0.09 ± 0.01 . Several 
radio pulsars with WD companions are known with compara- 
ble period and with highly circular orbits - e.g. 1855H-09 (Ryba 
& Taylor 1991) - as expected following stable RLOF from a 
low-mass giant. Probably the least far-fetched explanation of 
the eccentricity in VI 379 Aql is the presence of a third body in 
a substantially inclined orbit; and this might also explain some 
loss of angular momentum. 

5. CONCLUSIONS 

Case A RLOF, even when restricted to the classical 'conser- 
vative' model, shows a rich variety of behaviour, which we feel 
is often not emphasised enough. We identify 9 sub-classes, de- 
pending partly on whether the system evolves into contact (in 5 
different ways) or reaches reverse RLOF (in 4 different ways). 



Further subdivision depends on the evolutionary states reached 
when contact or reverse RLOF occurs. We can expect even 
more subclasses when non-conservative processes are modeled, 
as will certainly be necessary for extremes of high-temperature 
and low-temperature systems. 

For all but one of our selection of observed hot Algols we 
find an acceptable x^ when fitting the observed parameters to 
our library of conservative Case A binary tracks. It is encourag- 
ing to note that the worst outlier (AF Gem) lies near the lower 
boundary of the temperature range in which we expect the con- 
servative assumption to hold. The next largest x^'s come from 
two binaries with known third bodies, which may act to remove 
angular momentum from the inner orbit. 

Our selection of cool Algols shows significantly worse agree- 
ment between the observed systems and the conservative theo- 
retical tracks, suggesting the need for more free parameters in 
the modelling, such as mass and angular momentum loss. 

This data set of conservative Case A tracks has uses beyond 
an indiviual comparison of observed systems. With an esti- 
mate of the initial mass function and period distribution of bi- 
naries, it may be useful for population synthesis studies or for 
creating close binary-inclusive isochrones for stellar population 
studies. We hope to make these tracks available in early 2001 
on the Institute of Geophysics and Planetary Physics web site 
http://www.llnl.gov/urp/IGPP. 



This work was undertaken as part of the DJEHUTY project, 
which is developing a 3-D code to deal with hydrodynamical 
processes within stars, both single and binary. Most of the end- 
points to which our Case A systems evolved can be expected 
to have behaviour on a hydrodynamic timescale, and we hope 
to investigate them further in the future. We are grateful to our 
DJEHUTY colleagues for helpful discussions, and in particular 
to Don Dossa for helping with the parallelisation of the code. 
Work performed at LLNL is supported by the DOE under con- 
tract W7405-ENG-48. CAN is supported in part by a NPSC 
Graduate Fellowship. 
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Fig. 1 . — Case AD - Dynamic Contact. In the top panel of this figure, and in Figures 2 - 8, we show the evolution of both stars in an HR diagram. The track 
for *1 is always the initially more luminous. *2 the initially less luminous. Solid lines indicate periods where the binary is separated, dashed lines indicate periods 
where *1 is transfening mass to *2. We mark any transitions to the Hertzsprung Gap as H, and transitions to the Giant Branch as G. The bottom panel shows the 
mass transfer rate in logarithmic units of M/yr for the period in time during which mass transfer occurs. The initial parameters of this binary are logMio =0.15, 
qo = 0.50 and logPo/i^AMS = 0.15. The mass transfer rises rapidly to the dynamic timescale and the two stars come into contact. 
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Fig. 2. — Case AR - Rapid Contact. The initial parameters of this binary are logMio = 0.45, ^o = 0.15 and logPo/PzAMS = 0.10. The mass transfer rate rises 
rapidly to the thermal timescale and the two stars come into contact with a final mass ratio log g = 0. 11 , still well above unity. 
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Fig. 3. — Case AS - Slow Contact. The initial parameters of this binary are logMio = 0.45, log.jo = 005 and logPo/PzAMS = 0.20. Contact is avoided during the 
period of thermal scale mass transfer and a long period of nuclear timescale mass transfer follows. The run ends in contact with a final mass ratio, log q = -0.27. 
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Fig. 4. — Case AE - Early Overtaking. The initial parameters of this binary are logMio = 0.45, log^o = 0.05 and logPo/^ZAMS = 0.45. *2 gains so much mass 
that its evolution overtakes that of * 1 and *2 reaches the HG first. The run ends in contact with a final mass ratio of log^ = —0.23. 
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Fig. 5. — Case AG - Contact on Giant Branch. The initial parameters of this binaiy are logMio = 0.50, log^o = 0.10 and logPo/PzAMS = 0.55. The stars come 
into contact at a mass ratio of \ogq = -0.83. There is a brief period of separation at ? ~ 384 Myr with *1 in the HG. 
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Fig. 6. — Case AL - Late Overtaking. The initial parameters of this binary are logMm = 0.75, logi^o = 0.05 and logPo/^ZAMS = 0.60. In this run *1 loses so 
much mass in two periods of mass transfer that it eventually shrinks inside its RL and becomes a low mass helium burning core, logMi = -0.07. The run ends as *2 
crosses the HG and fills its RL at a very low mass ratio, log Cj = —\ .07. The brief period of separation between the two bursts of mass transfer is a feature common to 
all of our higher mass case AL/AN binaries. This feature occurs as * 1 exhausts hydrogen in the core, convection in the core shuts off and * 1 shrinks slightly inside 
its RL. At this point * 1 behaves as a "normal" massive terminal MS star, executing the classic hook at the end of the MS . When hydrogen is completely exhausted in 
the core and hydrogen shell burning begins *1 starts to cross the HG. During this rapid phase of envelope expansion, *1 quickly fills its RL again and mass transfer 
begins again. 
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Fig. 7. — Case AB - (Sub-type of AL). The initial parameters of tliis binaiy are logMio = 0.90, logi^o = 0.10 and logPo/PzAMS = 0.65. After igniting helium 
towards the end of the second burst of mass transfer, *1 shrinks inside its RL for awhile, but eventually re-expands and undergoes a third period of mass transfer. 
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Fig. 8. — Case AN - No Overtaking. The initial parameters of this binary are logMjo = 1.15, log^o = 0.05 and logPo/^ZAMS = 0.45. After two periods of mass 
transfer, * 1 becomes a helium star of mass log M[ =0.48. As *2 evolves towards the terminal MS, * 1 ignites carbon in the core. The ignition of carbon suggests an 
imminent supernovae explosion, and we conclude that *1 will supemovae before *2 fills its RL. 
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Fig. 9. — The classification of eacli of our individual binary nms into cases AD, AR, AS, AE, AG, AL, AB, AN, BD, BR, BL, and BN, as described in §3. Each 
square represents a slice through the data cube at constant Mjo- The upper left hand block is a slice at logMio = -0.05. The squares increase in the order one reads 
the page of a book, increasing in units of 0.05 and ending in the bottom right hand comer with logMio = 1.70. We do not show our results for logMio = -0.10 as 
very few of these binaries reached RLOF at an age younger then 20 Gyr. Within each square the x-axis represents increasing mass ratio in logarithmic units of 0.05 
from logij = 0.05 to 0.50. The y-axis represents increasing period in units of 0.05 from logPo/PzAMS = 0.05 to 0.75, where PzAMS is the critical period at which 
RLOF would occur on the ZAMS. The color of each dot represents the classification of the binary run according to the legend to the right of the plot. A white dot 
represents a Case A binary for which we could not determine a sub-class, a white dot outlined in black indicates a Case B binary for which we could not determine 
a sub-class. 
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Fig. 10. — The final mass ratio of each binary system. The ordinates are defined as in Figure 9. We define the final state to be the point where either a) contact is 
reached or seems imminent, b) reverse mass transfer begins, c) * 1 has ignited carbon or d) the binary is in a detached phase and reverse mass transfer is imininent. 
Binaries for which the final state could not be determined appear as white squares. Higher mass white squares are generally due to binaries which break down during 
thermal timescale mass transfer The color-coded final mass ratios are in logarithmic units. The legend is to be interpreted as follows: binaries with iogq < -1.0 are 
shown as yellow squares, those with -1.0 < log^ < —0.8 as orange squares, and so forth. 
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Fig. 1 1. — The fraction of a given binary's lifetime time .spent transferring mass from *1 to *2, /rlof- Tlie ordinates are defined as in Figure 9. The blanJi squares 
are those binaries which failed to converge for more then a few timesteps. The color-coded fractions are in logarithmic units. The legend is to be interpreted as 
follows: binaries with log/RLOF < ~3.5 are shown as yellow squares, those with —3.5 < /rlof < ~3.0 as orange squares, and so forth. 
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Fig. 12. — The residuals of observed binary fits. The hot Algols of Table 2 are shown in red, the cool Algols of Table 3 are shown in blue. The residuals of *1 are 
shown as circles, the residuals of *2 as squares. The residuals are defined in the sense A = fit— obs . 
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Table 1 
Summary of Binary Classification Scheme 



Case Defining Equations End State*! End State *2 End Geometry 



AD 


M > M/tiyn 


M 


M 


Contact 


AR 


M > M/fKH, 'contact -fRLOF(*l) < 0.10 X (ms(*1) 


M 


M 


Contact 


AS 


'contact -fRLOF(*l) > 0.10 X tMs(*l) 


m 


m 


Contact 


AE 


/h(*2)</h(*1) 


m 


H 


Contact 


AG 


'h(*1)<'h(*2) 


G 


H,G 


Contact 


AL 


fRL0F(*2)<'ctar,(*l) 


R,C 


H,G 


RLOF *2 


AN 


tebum(*l)<'RL0F(*2) 


SNe 


M,H,G 


Detached 



Table 2 
Hot Algols: observed, best fit and ZAMS parameters 



log log log log log log log log log 

Star P Ml <7m, q u,, Ti aj, Ti ctt- «i cts R2 tr U o-l, ^2 Ci Age 



0.030 

UCrB* 0.538 0.164 0.023 -0.420 0.022 3.767 0.015 4.170 0.009 0.694 0.008 0.436 0.011 1.430 0.060 2.510 0.130 

0.547 0.158 0.055 -0.481 0.055 3.761 0.049 4.187 0.031 0.703 ■■■ 0.433 0.032 1.403 ••• 2.567 

0.235 0.550 ■■■ 0.200 ••• 4.137 ••• 4.003 ■■■ 0.341 ■■■ 0.227 ••• 2.185 ••• 1.417 

AFGem' 0.095 0.063 0.015 -0.466 0.010 3.767 0.010 4.000 0.020 0.365 0.007 0.417 0.010 0.750 0.030 1.780 0.090 

0.050 

DM Per" 0.436 0.316 0.012 -0.547 0.011 3.920 0.010 4.260 0.020 0.677 0.008 0.653 0.009 2.000 0.060 3.280 0.100 

0.488 0.313 0.052 -0.474 0.051 3.908 0.039 4.248 0.034 0.715 ■■■ 0.650 0.031 2.016 ••• 3.243 

0.187 0.700 ■■■ 0.200 ••• 4.229 ••• 4.106 ■■■ 0.432 ■■■ 0.311 ••• 2.735 ••• 1.997 

VPup** 0.163 0.954 0.046 -0.277 0.068 4.360 0.060 4.420 0.040 0.724 0.024 0.799 0.020 3.850 0.250 4.200 0.160 

ATau^^ 



0.124 


0.732 


0.023 


-0.175 


0.010 


4.255 


0.020 


4.395 


0.020 


0.623 


0.010 


0.591 


0.011 


3.210 


0.030 


3.710 


0.149 


0.769 


0.055 


-0.201 


0.051 


4.242 


0.035 


4.384 


0.033 


0.650 




0.615 


0.031 


3.220 




3.720 


0.119 


0.950 




0.150 




4.369 




4.287 




0.581 




0.493 




3.593 




3.088 


0.538 


0.164 


0.023 


-0.420 


0.022 


3.767 


0.015 


4.170 


0.009 


0.694 


0.008 


0.436 


0.011 


1.430 


0.060 


2.510 


0.547 


0.158 


0.055 


-0.481 


0.055 


3.761 


0.049 


4.187 


0.031 


0.703 




0.433 


0.032 


1.403 




2.567 


0.235 


0.550 




0.200 




4.137 




4.003 




0.341 




0.227 




2.185 




1.417 


0.095 


0.063 


0.015 


-0.466 


0.010 


3.767 


0.010 


4.000 


0.020 


0.365 


0.007 


0.417 


0.010 


0.750 


0.030 


1.780 


0.124 


0.127 


0.052 


-0.313 


0.051 


3.776 


0.027 


4.023 


0.037 


0.414 




0.426 


0.032 


0.883 




1.897 


0.026 


0.400 




0.200 




4.038 




3.878 




0.254 




0.169 




1.612 




0.803 


0.312 


0.462 


0.029 


-0.409 


0.022 


4.064 


0.020 


4.300 


0.020 


0.643 


0.029 


0.763 


0.022 


2.490 


0.060 


3.680 


0.330 


0.497 


0.058 


-0.386 


0.055 


4.054 


0.039 


4.286 


0.033 


0.673 




0.757 


0.037 


2.516 




3.612 


0.120 


0.800 




0.150 




4.287 




4.200 




0.493 




0.402 




3.088 




2.554 


0.436 


0.316 


0.012 


-0.547 


0.011 


3.920 


0.010 


4.260 


0.020 


0.677 


0.008 


0.653 


0.009 


2.000 


0.060 


3.280 


0.488 


0.313 


0.052 


-0.474 


0.051 


3.908 


0.039 


4.248 


0.034 


0.715 




0.650 


0.031 


2.016 




3.243 


0.187 


0.700 




0.200 




4.229 




4.106 




0.432 




0.311 




2.735 




1.997 


0.163 


0.954 


0.046 


-0.277 


0.068 


4.360 


0.060 


4.420 


0.040 


0.724 


0.024 


0.799 


0.020 


3.850 


0.250 


4.200 


0.185 


0.927 


0.068 


-0.223 


0.085 


4.345 


0.065 


4.451 


0.045 


0.724 




0.793 


0.034 


3.784 




4.344 


0.117 


1.100 




0.100 




4.444 




4.395 




0.668 




0.610 




4.065 




3.755 


0.597 


0.276 


0.009 


-0.578 


0.007 


3.920 


0.030 


4.280 


0.040 


0.724 


0.016 


0.806 


0.007 


2.110 




3.690 


0.641 


0.321 


0.051 


-0.529 


0.050 


3.922 


0.048 


4.249 


0.048 


0.819 




0.803 


0.030 


2.281 




3.557 


0.254 


0.750 




0.200 




4.259 




4.138 




0.463 




0.341 




2.913 




2.185 


0.391 


0.362 


0.018 


-0.367 


0.039 


3.955 


0.020 


4.255 


0.040 


0.653 


0.019 


0.672 


0.018 


2.070 


0.060 


3.300 


0.387 


0.375 


0.053 


-0.417 


0.063 


3.949 


0.041 


4.245 


0.049 


0.670 




0.668 


0.035 


2.088 




3.268 


0.137 


0.700 




0.150 




4.229 




4.138 




0.432 




0.341 




2.735 




2.185 



16 



218 1.604 



635 11.394 



64 0.949 



114 3.379 



1.311 



99 2.918 



ZVul'" 0.391 0.362 0.018 -0.367 0.039 3.955 0.020 4.255 0.040 0.653 0.019 0.672 0.018 2.070 0.060 3.300 0.160 O 

107 0.776 ^ 
cc 

> 

Cd 

5' 
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Table 3 
Cool Algols: observed, best fit and ZAMS parameters 



Star 



log 
P 



log 

Ml 



log 



log 



log 



log 



R2 



log 



log 

Z.2 



Age 



S Cnc'^" 


0.977 


-0.638 


0.036 


-1.045 


0.001 


3.665 


0.010 


3.990 


0.010 


0.720 


0.004 


0.332 


0.004 


1.050 


0.050 


1.580 







1.033 


-0.605 


0.062 


-0.897 


0.050 


3.678 


0.039 


3.920 


0.036 


0.773 




0.326 


0.028 


1.210 




1.283 






-0.055 


0.150 




0.250 




3.832 




3.680 




0.148 




-0.144 




0.580 




-0.614 




RCMa'2" 


0.055 


-0.775 


0.049 


-0.801 


0.029 


3.630 


0.030 


3.860 


0.030 


0.025 


0.028 


0.196 


0.027 


-0.410 


0.160 


0.760 







0.152 


-0.580 


0.070 


-0.653 


0.058 


3.676 


0.048 


3.782 


0.039 


0.192 




0.214 


0.077 


0.041 




0.508 






-0.319 


0.000 




0.350 




3.751 




3.569 




-0.050 




-0.386 




-0.143 




-1.544 




RZ Cas^ 


0.077 


-0.137 


0.012 


-0.480 


0.010 


3.672 


0.020 


3.934 


0.005 


0.288 


0.007 


0.223 


0.008 


0.160 


0.080 


1.120 







0.109 


-0.192 


0.051 


-0.412 


0.051 


3.674 


0.043 


3.884 


0.036 


0.296 




0.233 


0.028 


0.239 




0.954 






-0.105 


0.150 




0.200 




3.832 




3.718 




0.148 




-0.102 




0.580 




-0.379 




TV Cas" 


0.258 


0.185 


0.014 


-0.393 


0.008 


3.720 


0.040 


4.020 


0.020 


0.517 


0.007 


0.498 


0.008 


0.860 


0.170 


2.030 







0.265 


0.134 


0.052 


-0.376 


0.051 


3.766 


0.060 


4.061 


0.037 


0.509 




0.503 


0.032 


1.033 




2.202 






0.097 


0.450 




0.200 




4.072 




3.924 




0.282 




0.183 




1.806 




1.014 




ASEri"" 


0.426 


-0.682 


0.018 


-0.968 


0.012 


3.720 


0.030 


3.930 


0.030 


0.340 


0.023 


0.196 


0.016 


0.470 




1.060 






0.542 


-0.587 


0.053 


-0.791 


0.051 


3.676 


0.048 


3.880 


0.048 


0.451 




0.186 


0.029 


0.558 




0.848 






-0.005 


0.150 




0.500 




3.832 




3.569 




0.148 




-0.386 




0.580 




-1.544 




TTHya''^ 


0.842 


-0.229 


0.132 


-0.646 


0.002 


3.680 


0.010 


3.990 


0.010 


0.769 


0.029 


0.290 


0.028 


1.200 


0.080 


1.500 







0.846 


-0.270 


0.141 


-0.633 


0.050 


3.652 


0.039 


3.999 


0.036 


0.759 




0.290 


0.040 


1.080 




1.529 






0.226 


0.200 




0.100 




3.879 




3.802 




0.167 




0.090 




0.803 




0.339 




AT Peg^' 


0.059 


0.021 


0.012 


-0.325 


0.008 


3.690 


0.017 


3.920 


0.005 


0.332 


0.006 


0.270 


0.006 


0.380 


0.070 


1.190 







0.090 


-0.025 


0.051 


-0.289 


0.051 


3.712 


0.031 


3.913 


0.036 


0.341 




0.266 


0.027 


0.483 




1.139 






0.054 


0.250 




0.250 




3.924 




3.750 




0.182 




-0.050 




1.014 




-0.149 




/3Per^" 


0.457 


-0.092 


0.026 


-0.663 


0.010 


3.650 


0.028 


4.100 


0.017 


0.544 


0.012 


0.462 


0.006 


0.630 




2.190 






0.554 


-0.085 


0.056 


-0.537 


0.051 


3.692 


0.047 


4.018 


0.036 


0.626 




0.466 


0.031 


0.973 




1.959 






0.154 


0.350 




0.200 




4.002 




3.831 




0.227 




0.150 




1.416 




0.578 




HU Tau** 


0.313 


0.057 


0.011 


-0.592 


0.008 


3.738 


0.012 


4.080 


0.034 


0.507 


0.004 


0.410 


0.005 


0.920 


0.050 


2.090 







0.414 


0.095 


0.051 


-0.405 


0.051 


3.733 


0.028 


4.078 


0.045 


0.594 




0.419 


0.031 


1.075 




2.102 






0.247 


0.450 




0.250 




4.072 




3.878 




0.282 




0.169 




1.806 




0.803 




TX Uma' 


0.486 


0.072 


0.014 


-0.606 


0.010 


3.740 


0.016 


4.110 


0.010 


0.627 


0.007 


0.451 


0.006 


1.170 


0.065 


2.300 







0.525 


0.096 


0.052 


-0.471 


0.051 


3.735 


0.030 


4.123 


0.031 


0.668 




0.461 


0.031 


1.229 




2.368 






0.266 


0.500 




0.250 




4.105 




3.924 




0.311 




0.183 




1.997 




1.015 




V1379Aql^ 


1.315 


-0.517 


0.021 


-0.873 


0.005 


4.490 


0.020 


3.650 


0.030 


-1.284 


0.076 


0.955 


0.036 


0.380 


0.110 


1.480 







1.372 


-0.577 


0.054 


-0.998 


0.050 


4.501 


0.043 


3.704 


0.047 


-0.997 




0.967 


0.044 


0.963 




1.702 






0.004 


0.250 




0.200 




3.924 




3.775 




0.182 




0.016 




1.014 




0.085 





,045 



180 



020 



090 



070 



030 



150 



040 



070 



4182 14.447 



19470 22.871 



3971 



502 



2078 



373 



5.382 



2.925 



3474 22.549 



2262 0.718 



2.252 



1135 15.972 



515 18.311 



8.113 



2585 10.138 



References. — (1) Cester et al. (1978); (2) Van Hamme & Wilson (1993); (3) Jeifery & Simon (1997); (4) Khalesseh & Hill (1992); (5) Maxted, Hill & Hilditch (1994a); (6) Maxted, Hill, 
& Hilditch (1994b); (7) Maxted, Hill, & Hilditch (1995a); (8) Maxted, Hill, & Hilditch (1995b); (9) Maxted & Hilditch (1996); (10) Popper (1973); (11) Richards, Mochnaclsi, & Bolton 
(1988); (12) Sarma, Vivekananda Rao & Abhyankar (1996); (13) Tomkin (1985); (14) Van Hamme & Wilson (1993) 



